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A simple microscopic principle of 'Valence bond (nearest neighbor singlet) amplitude maximiza- 
tion '(VBAM) is shown to be present in undoped and optimally doped cuprates and unify the very 
different orderings such as antiferromagnetism in the Mott insulator and the robust superconduc- 
tivity accompanied by an enhanced charge and stripe correlations in the optimally doped cuprates. 
VBAM is nearly synonymous with the energy minimization principle, ft is implicit in the RVB the- 
ory and thereby makes the predictions of RVB mean field theory of superconductivity qualitatively 
correct. 



The qualitative predictions of the RVB mean field the- 
ory of high To superconductivity in cuprates, including 
the symmetry of the order parameter has turned 

out to be in good agreement with experiments. In spite 
of its approximate character it has definitely put us in 
the correct Hilbert space by focusing on the key singlet 
correlations. On the other hand various orders, that at 
a superficial level appear to be unusual from the RVB 
theory point of view, have emerged both experimentally 
and theoretically. The aim of the present letter is to 
identify an unique cause behind the AFM order in the 
Mott insulator, the enhanced charge and spin stripe cor- 
relations and low energy spin fluctuations in the opti- 
mally doped superconducting cuprates 1^-^. A simple 
principle of 'valence bond (nearest neighbor singlet) am- 
plitude maximization' (VBAM) is shown to emerge as 
a fairly simple cause. This maximization respects cer- 
tain geometrical and dynamic constraint provided by our 
strongly correlated system. In the case of the Mott insu- 
lator, while super exchange interaction respects VBAM, 
an irreducible amount of local triplet fluctuations is gen- 
erated in order to satisfy certain geometric constraints. 
These triplet fluctuations interact and condense at mo- 
mentum (tt, ..tt) resulting in long range AFM order for 
any dimension > 2. In the optimally doped Mott insula- 
tor the constraints are modified by the charge dynamics. 
And VBAM continues as an approximate principle and 
produces new orders such as the robust superconductiv- 
ity and the enhanced quasi static charge and spin stripe 
correlations. We also argue that VBAM is responsible 
for the qualitative success of the RVB mean field theory. 

At a more fundamental level the 2 dimensionality, the 
one band character, strong correlation and closeness to 
a Mott insulating state are believed to be responsible 
for various interesting low energy electronic phenomenon 
in cuprates The VBAM hypothesis is an emergent 
consequence of the above and perhaps an useful guiding 
principle that at least helps us to separate cause from 
effects. One of the messages of this letter is that it is 
VBAM, rather than the development of charge or spin 
stripe correlations that is providing a mechanism for su- 



perconductivity: on the contrary whenever stripes are 
nurtured it is at the expense of superconductivity [^,0 . 

Our VBAM hypothesis is proved for the case of un- 
doped Mott insulator by simply using the principle of en- 
ergy minimization in the ground state. For the optimally 
doped Mott insulator various arguments are provided in 
support of this principle. The existence of the sharp 41 
meV triplet resonance [pT| at (tt, tt) at optimally doped 
bilayer YBCO and its counterpart at various dopings all 
the way up to zero doping is presented as an experimen- 
tal support for our VBAM hypothesis in the optimally 
doped regime. 

The phase coherent resonance of singlet bonds is at 
the heart of the RVB theory of superconductivity. An- 
derson's RVB wave function 

\RVB)^PG{Y.a.,,b\^)'^\0) (1) 

by construction has enhanced VB amplitude. It is a re- 
markably universal wave function in the sense, it can de- 
scribe an antiferromagnetically ordered Mott insulating 
state and the robust superconducting state of the opti- 
mally doped Mott insulator with equal ease for appropri- 
ate choices of the pair function . We believe that this 
universal feature is a consequence of the VBAM princi- 
ple. Also without changing the RVB wave function and 
consequently the superconducting property in a funda- 
mental fashion, the charge stripe and spin stripe order 
can be incorporated by a modulation of the pair function 
aij. 

To begin with we develop a heuristic picture of VBAM. 
In a free fermi gas, any short range singlet correlation 
contained in the ground state is a consequence of Pauli 
principle rather than interactions. In the large U Hub- 
bard model, accepted as a good model for our narrow 
band conducting cuprates close to half filling, every ele- 
mentary collision between two electrons tries to establish 
a nearest neighbor singlet correlations (a valence bond). 
That is, the virtual transitions to doubly occupied state 
on a given copper site lowers the energy (compared to the 
U = oo case) by the super exchange energy J Ri ^ and 
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stabilizes the spin singlet state rather than a triplet state. 
Elementary collisions, in addition to the Pauli principle 
induce singlet correlations. Close to half filling elemen- 
tary two body collisions are more frequent than free hop- 
ping of charges. The on site collision induced valence 
bond proliferation is at the heart of the VBAM principle. 

Let us try to understand the development of long range 
AFM order in the spin half Mott insulators in the light 
of the principle of VBAM. In an isolated pair of neigh- 
boring orbitals, the electron pair has a spin singlet (non- 
magnetic) ground state, through the super exchange in- 
teraction. An important question is how this local non 
magnetic singlet tendency manifests itself in a 2 or 3 di- 
mensional lattice. The spin half Heisenberg Hamiltonian 
with nearest neighbor interaction 

H = JJ2iS^■S,-^), (2) 

in terms of the underlying electron variables c's and the 
bond singlet operator blj = :^(CiT^]i ~ ''U'^iT-' ^^^^^ 
form 

H=^jY,bl,b^j, with n,|+n,j,^2 (3) 
(y) 

We have used the important identity [Q 

(S, . S, - in,n,) ^ bl^b,^ (4) 

The singlet number operator bljb^y has an eigen value of 
or 1 in our single occupancy subspace. In view of the 
above identity, in a translationally invariant ground state, 
what is maximized consistent with the lattice structure, 
is the valence bond amplitude. Thus in the nearest neigh- 
bor Heisenberg models minimization of the ground state 
energy is synonymous with maximization of the strength 
of nearest neighbor singlet bonds. This proves our VBAM 
hypothesis for the spin half Mott insulator. 

Even though diagonal in terms of the number opera- 
tors bljb^j, equation (2) is not really diagonalized, as the 
number operators themselves do not commute whenever 
one of the sites coincide: 

xSfc).S, (5) 

This non-commutativity propagates an irreducible mini- 
mum of bond triplet fluctuations in the ground state by 
making the ground state average {bljb^j)c < 1. In fact, 

{bljKj)G « 0.6391,0.5846 and = 0.5 respectively for Id 
chain, 2d square lattice and infinite d hyper cubic lattice. 
In the Id chain the bond triplet fluctuation is finite and 
manages to produce an algebraic AFM order. In 2d it 
increases further and is believed to produce a true long 
range order and in infinite d it increase even further to 
produce a perfect Neel order. The chiral operator ap- 
pearing on the right hand side of equation (5) also tells 



us that chiral fluctuations are also induced in the ground 
state. 

A long range magnetic order, when it occurs in the 
ground state, is an inevitable consequence of VBAM in 
the presence of the constraints provided by the lattice 
structure and the above commutation relation. Thus in 
a hyper cubic lattice for d> 2 

VBAM -I- geometrical constraints =^ AFM Order 

The nature of magnetic order is strongly lattice depen- 
dent. For the non bipartite 2d triangular lattice, what 
maximizes the bond singlet amplitude is a 120° structure 
with zero point fluctuations. The case of P doped Si in 
the insulating state is described by a 3d random lattice 
Heisenberg model. The nature of lattice constraint being 
very different long range magnetic order is believed to be 
absent resulting in a kind of singlet bond glass state. 

The role of the long range AFM order should not be 
also overemphasized in this context. As shown by Liang, 
Docout and Anderson ||l^ , the energy difference between 
a disordered spin liquid state with only short range AFM 
correlations and the the best variational state with long 
range AFM order is as small as one or two percent of the 
total energy: 

(Mam - jblA^h ^ 

At the level of variational wave function a small change in 
the long distance behavior of the pair function Oij takes 
us between a disordered and ordered ground state. Hsu 
p3[ also shows that the AFM order in 2d is a spinon 
density wave in a robust spin liquid state. Thus the de- 
velopment of long range AFM order in 2d is a result of 
a small final adjustment of the VB amplitude in a spin 
liquid state in the maximization procedure. 

We wish to argue that VBAM principle continues to be 
approximately valid for the optimally doped (S w 0.15) 
regime and produce, in the presence of the constraints 
modified by the hole dynamics, the robust superconduct- 
ing order and also the enhanced charge and spin stripe 
correlations. The discussion will be heuristic and also 
uses known results from the RVB theory. We will not 
consider the under doped regime, as disorder and long 
range coulomb interaction play important roles and sup- 
press superconductivity strongly and produce mesoscopic 
phase separation complication. 

The t-J model 

Hu = -t^{cl^Cj^ + h.c.) - J ^(S, • - ^n.nj), 

takes the form 

Htj = ^(cLc,, + h.c.) -JJ2 bl,b.,^ (7) 
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with the usual constraint n^i + riii ¥^ 2. We would like 
to see if VBAM continues to be valid in the presence of 
the single electron hopping term in the optimally doped 
case. 

From energy minimization point of view the double oc- 
cupancy constraint limits the single particle kinetic en- 
ergy gain per site to Ake ~ tzd compared to the larger 
super exchange energy gain Ase ~ Jz(l — 5^). For 
cuprates, with a co-ordination number z = 4 and j « 2, 
when S ^ 0.15 the ratio w 3. From the energy 

considerations given above, the VBAM principle contin- 
ues to be important in the optimally doped regime. The 
above rough estimate is in agreement with more accurate 
estimates using variational monte carlo studies. 

The nature of constraints on valence bond proliferation 
has changed now; we will see how it can stabilize new low 
temperature phases such as the robust superconductiv- 
ity along with quasi static charge and spin stripe correla- 
tions. Maximizing valence bond amplitude is an impor- 
tant step towards establishing a robust superconducting 
state. The next important step is the development of 'in 
phase resonance' of the valence bonds (or the zero mo- 
mentum condensation of the charged valence bonds) in 
the ground state. This is what is precisely achieved in the 
RVB mean field theory. In this sense the VBAM princi- 
ple is satisfied by the RVB mean field theory. While the 
original RVB mean field theory emphasized the extended- 
s mean field solution, Kotliar's identification of d^2_y2- 
wave RVB solution as a lower energy mean field solution 
made the RVB theory closer to experiments in terms of 
the symmetry of the order parameter. 

After the original RVB mean field theory, and some 
experimental developments of that time, the inclusion 
of the on site constraint was suspected to produce a 
large phase fiuctuations and in the even remove the fi- 
nite temperature K-T transition in an isolated C'u02 
layer. However, recent experiments have strongly 
supported the one layer d-wave superconductivity with a 
large 95K; in addition the RVB Ginzburg Landau 
functional derived by Anderson and the present author 
1^,^, in the RVB gauge theory approach did not show 
any singular effect on the GL coefficients due to the on 
site constraints. Very recently Lee 0, using RVB gauge 
theory has provided a rather convincing non perturbative 
analysis supporting the local stability of the d-wave solu- 
tion. 

What the experiments and also various theoretical 
clues have been telling us is that the effect of the strong 
correlation (on site constraint) does affect superconduc- 
tivity in unexpected fashions by encouraging quasi static 
charge stripe and spin stripe order. Within the VBAM 
principle we can understand it in the following heuris- 
tic fashion: the super exchange term tries to segregate 
the holes so that they have more fiuctuating Mott insu- 
lating region where super exchange gain (or equivalently 



the bond singlet amplitude) can be maximized. If the 
hole segregation has the form of a Id charge stripes, in 
addition to VBAM, kinetic energy gain also results from 
coherent charge delocalization along the charged stripes 
for the following reason. In view of the Brinkman-Rice 
phenomenon and proliferation of non self retracing paths 
in 2d, 'holes' can not maximize the kinetic energy gain 
by coherent delocalization in 2d. The quasi Id fluctuat- 
ing stripes on the other hand provides maximum of self 
retracing paths (along the stripes) for the holes thereby 
gaining coherent charge delocalization energy. 

Since the quasi static charge stripe formation also cre- 
ates coherently fluctuating Mott insulator region, VBAM 
leads to enhanced quasi static antiferromagnetic correla- 
tions (incommensurate order). Thus charge stripe and 
spin stripes are intimately related. In the optimally 
doped region, the development of charge and spin stripe 
correlations arise from finer adjustments of a robust su- 
perconducting state in order to satisfy the VBAM princi- 
ple. These adjustments can be already done at the level 
of RVB mean field theory in the following fashion both 
in the normal state and superconducting state: a) devel- 
opment of spontaneous anisotropic valence bond ampli- 
tudes jAj;! ^ I A J, I and b) one dimensional spatial 
modulation of |Ay | either along the x or y axis Jl8| , p^ . 

The first case corresponds to stripes having a orienta- 
tional order (along the a or b axis) and no spatial order. 
In the normal state this will be a nematic metal with 
Pa 7^ Pb and other wise a non fermi liquid state. This 
state will have an enhanced magnetic correlations around 
(tt, tt), in view of the spin localization caused by the quasi 
static stripe formation, This anisotropy can continue into 
the low temperature superconducting phase and will re- 
sult in a real mixing of d^2_y2 and extended-S state. Per- 
haps this phase is already seen j20) in the under doped 
LSCO - we think that the conducting charge stripe glass 
state is an anisotropic metallic state (nematic metal) . As 
explained elsewhere [|j coupling to octahedral rotations 
and displacements in LSCO further enhance the charge 
stripe correlation at the expense of superconductivity. 

The second case corresponds to a spatial ordering of 
the stripes, which can also start in the normal state, in- 
dependent of the low temperature superconducting state. 
Close to the magic filling of | this is believed to happen 
in doped LCD. 

In both cases there will be domain formation and the 
transition from the isotropic metallic phase to the above 
anisotropic phase can be a second order phase transition. 
In the optimally doped region, even in the most favorable 
case of LSCO, the stripe correlations are not stabilized 
into a true long range order. However, they remain as 
additional correlations in the ground state there by re- 
ducing superconductivity. 

At the level of wave function modifications the above 
secondary stripe orders appear as modifications of the 
pair function atj without changing the RVB supercon- 
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ducting wave function in a fundamental fashion. 

We will discuss briefly the question of low energy spin 
fluctuations in optimally doped cuprates. It is an ex- 
perimental fact that in the singlet dominated cuprates 
there are strong low energy spin fluctuation phenomenon 
in k-space particularly around the (tt, tt) region. Among 
them a remarkably neat and unique phenomenon occurs 
in the bilayer cuprates: in neutron scattering a resolution 
limited 41 meV peak corresponding to a triplet exci- 
tation with momentum centered around (tt, tt) is seen be- 
low the superconducting Tc- These excitations have very 
little dispersion over a wide momentum interval around 
(tt, tt), suggesting that we can create spin triplet wave 
packets of size comparable to the lattice parameter as 
nearly exact eigen excitations. That is, the triplet exci- 
tation is predominantly made of nearest neighbor (Eder 
pl| also has a nearest neighbor triplet bond picture for 
the triplet resonance) triplet bonds. Availability of large 
valence bond amplitude is a prerequisite for being able to 
create this excited state. It tells us that large amplitude 
valence bond exist in the ground and they can be indeed 
converted into triplets carrying a momentum of (tt, tt) 
and energy of about 41 meV. And 41 meV is the sta- 
bilization of bond singlet energy in the superconducting 
state relative to the normal state. Indeed various authors 
p2| , p3t have extracted the superconducting condensation 
energy from the spectral properties of the 41 meV peak. 
It should be also mentioned that interlayer pair tunnel- 
ing adds further stability to the valence bonds in bilayer 
cuprates there by supporting our VBAM principle. 

We suggest that the triplet resonance around (tt, tt) is 
what carries the memory of the Mott insulator. It con- 
tinues to persist as we go to lower dopings in YBCO and 
the energy of this peak decreases linearly with the cor- 
responding superconducting transition temperature. In 
a Mott insulator, where the superconducting Tc vanishes, 
this triplet excitation around (tt, tt) becomes soft (in view 
of the existing triplet condensate) and becomes a Gold- 
stone mode of the antiferromagnetic order. This also tells 
us that the maximization of the valence bond amplitude 
is easier in a doped Mott insulator than in an undoped 
one. The hole dynamics in some sense decreases the con- 
straints provided by the lattice, by effectively converting 
the valence bond operators into more of 'unconstrained' 
bosons. 

Thus the various antiferromagnetic fluctuations includ- 
ing the narrow resonance peak tells us about the growing 
antiferromagnetic correlation in the ground state, aris- 
ing from the stabilization of the bond singlets. It is in 
this sense the antiferromagnetic fluctuations are effects 
of a deep and growing bond singlet tendency, forced by 
some geometrical and dynamical constraints, rather than 
some thing that provide cooper pairing at low energies 
p4|. Spin fluctuations are effects rather than causes of 



the singlet dominated superconductivity phenomenon. 

Thus the maximization of the bond singlet amplitude 
can be thought of as a driving force behind the vari- 
ous anomalous correlations one sees in cuprates includ- 
ing the robust high Tc superconductivity. The emergent 
quasi static charge and spin stripe correlations and su- 
perconductivity mutually adjust themselves when some 
parameters such as doping or temperature are changed; 
they are not to be thought of driving each other. The 
fundamental driving force is VBAM which is implicit in 
the RVB theory. 



[1] 

[2] 
[3] 

[4] 

[5] 

[6] 
[7] 
[8] 



[9 

[lo; 

[11 
[12; 

[is: 

[14 

[is: 
[16: 

[17 

[is: 

[19 
[20 
[21 
[22' 
[23^ 
[24' 



1. J.G. Bednorz, and K.A. MuUer, Z. Phys. B 64, 189 

(1986) 

P.W. Anderson, Science 235 1196 (87) 
P.W. Anderson, The theory of high temperature super- 
conductivity (Princeton University Press, NY, 1996) 
G. Baskaran, Z. Zou and P.W. Anderson, Sol. State 
Commm. 63 973 (87) 

G. Kotliar, Phys. Rev. B 37 3664 (88); I. Affleck and J. 
B. Marston, Phys. Rev. B 37 3774 (88); H. Fukuyama, 
Prog. Theor. Phys. (Supplement) 108 287 (92) 
Zaanen, J and Gunnarsson, O., Phys. Rev. B 46 7391 - 
1394 (1989) 

V. Emery, S.A. Kivelson and O.Zachar, Phys. Rev. 56 
6120 (1997) 

J.M. Tranquada et al. Phys. Rev. Lett., 78 338 (1997); 
J.M. Tranquada, Proc. of Univ. of M iami Conf. on High- 



:ond-mat/990345S 



TcSuperconductivity (Jan, 1999), 
V.J. Emery, S.A. Kivelson and J.M. Tranquada, 
Natl. Acad. S ci. (USA). 96 8814 (99) 



Proc. 



ond-mat/991016l| (to appear in IJMP Let- 



G, Baskaran, 
ters B, 2000) 

S.R. White and D.J. Scalapino, Phys. Rev. B 60 R753 
(1999); Phys. Rev. Lett., 84 3021 (00) 

H. F. Fong et al., Phys. Rev. Let., 75 316 (95) 

S. Liang, B. Doucot and P.W. Anderson, Phys. Rev. 

Lett., 61 365 (1988) 

T. Hsu, Phys. Rev. B 41 11379 (90) 

K.A. Moler et al.. Science 279 1193 (1998); C. Panagopo- 

los et al., Phys. Rev. Lett. 79 2320 (1997) 

G. Baskaran and P.W. Anderson, Phys. Rev. B37 580 



M. Drzazga et al., Z. Phys., B74 67 (1989); 

D.H. Lee, Phys. Rev. Lett., 84 2694 (00) 

M. Vojta and S. Sachdev, Phys. Rev. Lett., 83 3916 (99); 

see also, C. Cas tellani et al., Phys. Rev. 75 4650 (1995). 

J.H. Han et al., |cond-mat/000604( 

A.W. Hunt et al., Phys. Rev. Lett., 82 4300 (1999); 

R. Eder, Phys. Rev. 59 13810 (99) 

D. J. Scalapino and S.R. White, Phys. Rev. 58 8222 (98) 

E. Demler and S.C. Zhang, Nature 396 733 (98) 
D. Pines, Phys. Rep. 250 329 (95) 



4 



